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Abstract 

We study perturbations of black iioles absorbing dark energy. Due to tlie accretion of dark energy, 
the black hole mass changes. We observe distinct perturbation behaviors for absorption of different 
forms of dark energy onto the black holes. This provides the possibility of extracting information 
whether dark energy lies above or below the cosmological constant boundary w = ~1. In particular, 
we find in the late time tail analysis that, differently from the other dark energy models, the accretion 
of phantom energy exhibits a growing mode in the perturbation tail. The instability behavior found 
in this work is consistent with the Big Rip scenario, in which all of the bound objects are torn apart 
with the presence of the phantom dark energy. 
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There has been growing observational evidence showing that our universe is accelerated expanding driven 
by a yet unknown dark energy (DE) [11,0,01. The leading interpretation of such dark energy is a cosmological 
constant with equation of state w — —1. More sophisticated models have been proposed to replace the 
cosmological constant by either relating the dark energy to a scalar field called quintessence with w > —1, 
or to an exotic field called phantom with w < —1. But it is doubtful that there is a clear winner in sight to 
explain the nature of dark energy at the moment. Recently, extensive analysis found that the current data 
favors dark energy models with equation of state in the vicinity of w = —1 straddling the cosmological 
constant boundary. This observational value was pinned down through large scale surveys from CMB, large 
scale structure and SNIa observations. 

At present the equation of state is the only antenna to learn the microscopic nature of dark energy. Its 
observational value is particularly important to determine whether the dark energy is of the phantom type, 
quintessence type or cosmological constant. Besides the available observational methods, it is of great interest 
to devise other complementary tools to measure the values of cj. In this work we use quasinormal modes 
(QNM) of a black hole to investigate the nature of DE equation of state. As is well known, if the dark energy 
is modeled as a background cosmic fluid, the flux of dark energy accreted by the black hole will change its 
mass. Now, it is reasonable to take into account these imprints in perturbations around a black hole and to 
extract information about w. An attempt in this direction was studied in 6| by exploiting the gravitational 
wave radiated from a binary of supermassive black holes. It is expected that the binaries observed with LIGO 
or LISA can distinguish, in the nearly future, more accurately values of w of dark energy. 

To extract information from the gravitational wave observation and pin down the value of w of dark energy, 
we need accurate waveforms on the perturbations around black holes. There has been great progress in 
studying perturbations around black holes recently. In asymptotically flat spacctimes it is possible to get a 
schematic picture regarding the dynamics of waves outside black holes. After the initial pulse, the perturbation 
will experience a quasinormal ringing, which describes the damped oscillations under perturbations in the 
surrounding geometry of a black hole with frequencies and damping times of oscillations entirely fixed by the 
black hole parameters. The quasinormal modes is believed as a unique fingerprint to directly identify the black 
hole existence. Detection of these QNM is expected to be realized through gravitational wave observation 
in the near future 7|. At late times, quasinormal oscillations are swamped by the relaxation process. This 



relaxation is the requirement of the black hole no hair theorem 



8|. In Anti-de Sitter and dc Sitter spacetimes. 



perturbations around black holes have been shown as a theoretical testing ground to get deeper understandings 



of the AdS/CFT 
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IJ] correspondences. More recently, it has been also argued 



that perturbations in black hole backgrounds is useful to extract information on black hole phase transitions 
In the universe filled with dark energy modeled as scalar field, the action has the form 



(1) 



where /i is the mass of the scalar field, and wc have taken the metric sign (— , +, +, +)^ If the dark energy is 
described by the phantom field, the kinetic term in the action has the positive sign ^18] , Varying the action 
with respect to ip, we can obtain the wave equation in the curved spacctime 
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(2) 



where the '+' sign describes the the phantom field while the ' sign describes the quintessence field. The 
wave equation of the quintessence field is the same as the massive scalar field discussed in 19 1 . 



Considering next the effect of the accretion of dark energy onto the black hole, the black hole mass changes 
at a rate [5| 

dM 



dt 



47rAA/>oo(l + u;), 
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where A ~ 5.6 is the numerical factor, which define the energy flux of DE onto the black hole, poo is the energy 
density of dark energy far away from the hole. The black hole is no longer static due to the absorption of dark 
energy, unless the dark energy being cosmological constant. Accreting the quintessence field with w > —1, 
the black hole mass increases. However the black hole mass decreases during the accretion of phantom energy 
with w < — 1. It was argued that the black holes are not torn apart, but disappear by the Big Rip due to 

n 

the accretion of phantom energy [5|| . Other discussions on the change rate of the black hole mass due to the 



absorption of dark energy can be found in [20| . What would be the fate of the perturbation around black holes 
absorbing dark energy? To answer this question we shall employ the formalism^first developed in Ref. 
Some other studies on the time-dependent background QNMs can be found in 



22| 



2ii. 
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To describe the time-dependent black hole background, wc start with the Vaidya metric 

ds^ = -fdv^ + 2cdvdr + {dO^ + sin^ edc/}'^), (4) 

where / = 1 — '^'^^^^'^ and M{v) is an arbitrary function of time. The coordinate v is usually called "advanced 
time" and c = 1 describes black hole with ingoing radial flow. The horizon of the Vaidya black hole is inferred 



from the null hypcrsurfacc condition g^'^ 



and f{r+,v) = 0. In this case, r+{v) satisfies 



r+ - 2M (v) - 2cr+r+ = 0, 



(5) 



where r+ = 

Similarly to the static approach, we introduce the tortoise coordinate as 

r* = r + -!-ln(r - r+), 
where k is the surface gravity. Since r» = r» (r, r_|_(w)) ,dv^dt + dr^,, one has 

— - (l 
dv \ dv J dt 
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Using the Vaidya metric and the tortoise coordinate, the wave equation ^ in the time-dependent black 
hole background can be written as 
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^- = 0, 
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where before //^ term the '+' sign describes the phantom field while the '— ' sign is for the quintessence 
field. Since the Vaidya metric is spherically symmetric, above equation thus is separable and the radial wave 
equation can be written as follows 



drl dv^dr^ 
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dr^. 
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where 
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(12) 



with -J ~ = 1 + 2K{r^-r+} ■ Here in the effective potential V, the '— ' sign is for the phantom field and the 
'+' sign is for the quintessence field. For the convenience of later numerical calculation, we adjust k in a way 
that Eq.® for each multipole moment becomes the standard wave equation near the horizon r_)-(t;o) 7|. The 
boundary condition at the horizon reads 



lim £2 = lim Si 
r ^ r±{vo) r ^ r±{vo) 

V Vq V ^ Vq 



lim F = 0, 

r r±{vo) 
V ^ vo 



(13) 
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which requires k = • 




FIG. 1: QNM behaviors of black holes absorbing dark energy. In plotting the figure we have taken parameters I = 2, 
Poo = 0.001, A ~ 5.6, ^^ = 0.001. 
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FIG. 2: The real and imaginary parts of quasinormal frequencies when black hole absorbs different types of dark energy. 
We have taken parameters I — 2, poo = 0.001, A ~ 5.6 and jj, — 0.001 in plotting the figure 



In order to simplify we make a variable transformation u = u(r^,v^) and v = v^, where the curve 
u(r^,v^) = constant is determined by the equation 

aw* 2c 



In addition, when c = 1 and £2 ^ 0, we have back the usual null coordinate u — > — 2r*, similar to the static 
Using the new variables u and u, the wave equation changes to 



case 
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In the static limit, £2 = £1 = 0, and w = w.^ — 2r*, thus the radial equation returns to its usual form 7| 

1/ 



Now for the Vaidya spacetime, u(r*,u*) satisfies 



du 



2c 



= 0. 







(16) 
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and Eq. p5)) simplifies to 



52 / 9u 



= 0. 



(18) 



Numerical solution of (|18p proceeds by first integrating w(r*,u*) according to (|14|) . and then using next 
the finite difference method developed in 



21| to solve the wave equation in the Vaidya metric. Denoting 



^ = *(-Ui,Wj), and 



i9u Au ' 9w Aw 

we can discrete ([T8| into 
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Taking n = = *s = *w = ^^-Lj+i, and using ^i^j = *ii:±*£, we have 
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du\ ^^E + '^w 
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FIG. 3: Late time evolution of the perturbation of massive scalar field in the stationary spacetime. We have taken 
w = —1, and different values of the field mass fj, — 0.001, 0.005 and 0.01 respectively. 
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FIG. 4: Late time evolution of the perturbation in around the black hole absorbing quintessence and phantom fields. 
We have taken different values of the field mass as ^ = O.OOf , 0.005 and O.Of respectively. 

Now we report our numerical results on the perturbations around the black hole with accretion of different 
types of dark energy. Fig.l illustrates the QNMs behavior. The solid line represents the case of accretion of 
phantom type dark energy with w = —1.1, the dashed line represents the dark energy being the cosmological 
constant, and the dotted line shows the result of absorbing quintessence field with w ~ —0.9. The details 
of quasinormal frequencies iOfi^LUj as functions of v are shown in Fig. 2. When the black hole absorbs the 
cosmological constant, its mass docs not change so that both the real and imaginary parts of quasinormal 
frequencies remain constant, in agreement with the static black hole cases. The black hole mass increases 
due to the accretion of quintessence field. Differently from the stationary black hole case, we observed that 
both the real part and the absolute value of the imaginary part of quasinormal frequencies decrease with the 
time evolution. However due to the accretion of phantom field, which causes the black hole mass to decrease, 
both the real part and the absolute value of the imaginary part of quasinormal frequencies increase with time. 
Comparing the values of the imaginary parts of QNM, we found that \ujj\ is smaller for black holes absorbing 
quintessence than phantom field. This explains what we observed in Fig.l, in which the effect of perturbation 
can last longer in the black hole background due to the accretion of quintessence type dark energy. 

We have also investigated the late time tail behavior of the perturbation in the black hole due to the 
accretion of dark energy. If the dark energy is cosmological constant, the result is shown in Fig. 3. Since 
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FIG. 5: Behaviors of the effective potential V{r) at the moment v = 100 for different values of fx for the accretion of 
phantom and quintessence respectively. In plotting the figure, we have taken I = 0, pcx, = 0.001, A ~ 5.6. 

the black hole mass do not change in this case, our result gives the objective picture of the late time tail 
behavior for the massive scalar field perturbation in the stationary black hole background. Differently from 
massless scalar field perturbation, we observed oscillations in the tail of massive scalar field perturbation. The 
frequency of the oscillation increases with the value of mass of scalar field fi, which confirms the argument 
proposed in [IS | . 



The late time tails of perturbations around black holes absorbing quintessence field and phantom field are 
shown in Fig. 4. The accretion of quintessence field has the similar oscillatory decay behavior in the late 
time tail as that observed for black hole absorbing dark energy with w = — 1. Interestingly, the accretion of 
phantom onto the black hole presents us completely different behavior in the late time tail. Instead of decay, 
we saw the growing modes in the t ail. The same phenomenon has also been observed for static black holes 
surrounded with phantom type DE 



24l | . Here we found that the growing appears earlier and faster when fi is 
bigger. The growing modes will make the black hole unstable. Thus from the tail of the perturbation we saw 
that it is not peaceful for the black hole to disappear by the Big Rip due to the accretion of phantom energy 
as argued in 5[. In the later moment when the black hole has accumulated enough phantom energy, it will 
"explode" . To find physical reason for this "explosion" , let us examine the behavior of the effective potential 
in the wave equation. Since the black hole mass will change due to the accretion of phantom or quintessence 
fields, the effective potential is time-dependent. Choosing v = 100, the behaviors of the potential are shown 
in Fig. 5. For the accretion of phantom, the effective potential outside the black hole approaches a negative 
value, which is different from the case with accretion of the quintessence. Actually this can be seen from the 
expression of the potential (jl2[) . whose behavior does not change with the time evolution. The existence of 
the negative part in the effective potential can form true bound state which leads to the growing modes 25 1. 

In summary we have investigated perturbations around black holes absorbing dark energy. It is well known 
that the mass of a black hole changes due to the accretion of dark energy of quintessence and phantom 



types 



5|. The study of perturbation in this black hole background is more complicated than in the usual 



2l| . we have observed different QNM behaviors for 



stationary situation. Using the formalism developed in 
the accretion of different types of dark energy onto the black hole. The QNM results discussed here are 
sufficient to illustrate the possibility to distinguish whether dark energy lies above or below the cosmological 
constant boundary w = — 1 in the future observation by exploiting the perturbations around black holes. 
However since the change rate of the mass due to the accretion of dark energy is extremely slow, to extract 
very accurate information of w from gravitational waves emitted from black holes while avoiding background 
noise still remains a challenge. 

The late time tail behavior of the perturbation around black hole absorbing dark energy presents us interest- 
ing results. Instead of the anticipated oscillatory decay in the background of black hole absorbing cosmological 
constant and quintessence, the accretion of phantom energy onto the black hole exhibits a growing mode in 
the perturbation. This growing tail tells us that due to the accretion of phantom energy the black hole does 
not disappear peacefully as argued in [5|. Instead, the black hole will explode after getting enough phantom 
energy. The result is consistent with the Big Rip scenario, i.e, all of the bound objects are to be torn apart 
in the presence of the phantom dark energy. 
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